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We calculate on the lattice the interface tension in the SU (2) pure gauge theory in d = 2 + 1 at high temperature. 
The result is compared to the perturbative prediction. The agreement confirms applicability of the perturbation 
theory in this case. 



1. Hot QCD and Zjv symmetry 

Studies of QCD at high temperature attracted 
considerable attention recently (see, e.g., [1,2] and 
refs. therein). One of the interesting questions 
concerns the existence and physical significance of 
the so-called Zjv (N = 3) phases of quark-gluon 
plasma at high temperature (T > Tc). These 
phases appear naturally in Euclidean formulation 
of QCD. In such a formulation the system lives 
in a 4-dimensional box. The extent of the three 
spatial dimensions is large compared to the phys- 
ical scale and the extent of the fourth (Euclidean 
time) dimension is l/T. The boundary condi- 
tions in this time direction are periodic (for gluon 
fields). 

Consider a hypothetical QCD without quarks 
and generalize it to an SU (N) gauge theory in 
d space-time dimensions. An order parameter 
distinguishing between the Zjv phases is the value 
of the Polyakov loop: a path ordered exponent 
along a path in the time direction at fixed spatial 
coordinate x: 

P(x) = j^TiPexY>i^igl^' rftAo(t,a!)|, (1) 

where Af, = A^^^T"- is the N x N matrix of the 
gauge potential. It is gauge invariant due to peri- 
odic boundary conditions provided that the gauge 
transformation U(x,t) is also periodic in t. How- 
ever, there are also gauge transformations which 
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are periodic up to an element of the center of the 
gauge group: 

U(x,t + l/T) = zU(x,t), (2) 

where z = I exp{i2Trk/N), A; = 0, 1, 2, 1; 
7 is a unit N x N matrix. The periodic b.c. are 
not affected by such a gauge transformation (if 
there are no fermions) as z commutes with all ma- 
trices in the gauge group. However, the Polyakov 
loop acquires a phase exp(2i7rA;/#). 

If one integrates out all degrees of freedom ex- 
cept for the Polyakov loops the resulting effective 
theory of the complex scalar field P(x) in d— 1 di- 
mensions will have a global Zjv symmetry: P 
exTp(i2'iTk/N)P . It turns out [3] that at high T this 
system is ferromagnetically ordered: (P) ^ (be- 
cause an effective ferromagnetic coupling grows 
with T). The value of (_P) can be related to the 
free energy _F of a static quark inserted in the hot 
gluon plasma: (_P) = exp(— _F/T). Thus confine- 
ment corresponds to (_P) = 0, i.e., F = cxd. In 
this way one can relate the breaking of the Zjv 
symmetry to the deconfinement transition. 

2. Interfaces and perturbation theory 

There are N possible complex directions where 
the vacuum expectation value (_P) can point. If 
two parts of space are occupied by domains with 
different values of (_P) an interface occurs between 
the domains. Such interfaces might have interest- 
ing cosmological consequencies [4] . 

Perturbation theory can be used to calculate 
properties of such interfaces at very high tem- 
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perature when the effective gauge coupiing be- 
comes smaii. Integrating out quadratic fluctu- 
ations around constant homogeneous configura- 
tion Aq one obtains the effective potential as a 
function of Aq. Consider, for example, an SU(2) 
gauge theory. One can rotate the constant field 
Aq to some direction in the matrix space, say, 
Aq = AqT3/2. Then the effective potential has 
the form as shown in Fig. 1 [5]. 
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Figure 1. Vefi- 



The periodicity of Veff refiects the Z2 sym- 
metry of the effective theory of the Polyakov 
loops: P —P. Indeed, for our constant Aq: 
P = cos (gAl/2T). At high T the Z2 symmetry 
will be spontaneously broken: minima of Veff cor- 
respond to _P = ±1. 

One can calculate the interface tension between 
two phases by considering a soliton-like solution 
which starts in one minimum, goes through the 
barrier and ends in another minimum. The ac- 
tion on such a trajectory would give the interface 
tension. In c? = 4 one obtains [6]: a = cT^/g, 
where c is a numerical constant. This result is 
easy to understand: the width of the interface is 
of order l/gT - Debye screening length, and the 
energy density inside the interface is of order T^. 

One has to take perturbative results at high 
temperature with caution, however. Corrections 
from higher orders of the perturbative expansion 
contain infrared divergencies. For example, it was 
suggested in [7] that infrared divergencies might 
result in vanishing interface tension. 

The aim of this paper is to study a pure gauge 
theory nonperturbatively (on the lattice), deter- 
mine the interface tension at very high tempera- 
tures (T ^ Tc) and compare these results with 
perturbation theory. 



3. Lattice study 

We study SU (2) pure gauge theory in 
c? = 2 -|- 1. The gauge variables — SU(2) matri- 
are defined on the links of a simple cubic 
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Ly X Lf lattice. The action is a sum 1 



the 



plaquette products Dp of the gauge matrices: 



^ = E(l-^Trnp) 
p 



(3) 



The partition function is Z = exp(— /JS'). 

To enforce an interface we use a twist [8]. The 
action of a system with the twist differs from 
the untwisted action in that the contribution of a 
set of plaquettes pierced by a line in y direction 
changes sign. One can check that this is equiv- 
alent to antiperiodic b.c. in x direction on the 
Polyakov loops. 

We simulate the system with and without the 
twist and measure the average action in the two 
cases. The difference AS is related to the excess 
free energy AF = — lii(Ztvi/Z) due to the inter- 
face and hence to the interface tension a: 



AS 



OAF 

dp 



d_ 

dp 



(tA 



(4) 



where A = Lya \s the "area" of the interface, 
T = l/(Lia) is the temperature in physical units 
and a is the lattice spacing. 

We calculated the interface tension for this 
model in continuum to the leading order in g^: 



ao = 5.104. 



(5) 



The interface width is of order g^/T. The g^ is 
the (dimensionful) bare coupling constant for the 
SU(2) gauge theory in continuum. It is related 
to the lattice parameter /3 as: /3 = A/^g'^a). The 
dimensionless small expansion parameter of the 
perturbation theory is g'^ jT , which in terms of 
the lattice parameters is: ALt/ (3. 

We can now substitute (5) into (4) and express 
A, T and g through the lattice parameters Ly, Lt 
and [3 to get the leading large [3 behavior of AS 
in perturbation theory: 
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Figure 2. The value of a = 4_L( ^S\/]3 / Ly as a 
function of at a given Lt and a linear extrap- 
olation to /3 = cxD. (For comparison: 1/ (3c ~ 0.29 
[9]). 

We perform measurements at different values 
of [3 at given Lt = 2,3,4. The spatial sizes L^, 
Ly are chosen after test runs to study finite size 
dependence. For Li = 2 and /3 = 15 we chose 
Lx = 36 and Ly = 12 and then scaled the sizes 
with the interface width (roughly as \/Ltj3). We 
plot AL^J"^ ^S\/]3 1 Ly as a function of 1//3 and ex- 
trapolate to /3 = cxD (see, e.g.. Fig. 2). The result- 
ing tto depends on Lt which is the cutoff 1/a in 
units of the temperature T . The continuum limit 
is Lt oo. It is possible, however, to calculate a 
to the leading order in g'^ /T ~ 1//3 in perturba- 
tion theory for finite Lt- We compare our results 
to perturbation theory at each Lt in Fig. 3. 

The extrapolation to 1//3 = (i.e., T = cxd) 
in Fig. 2 has an ambiguity as we do not know 
the form of the perturbative correction to a in 
(5). The linear 1//3 extrapolation we used as- 
sumes this correction to be 0(gT^^^). While we 
do not know its precise form, a simple analy- 
sis shows that this correction could contain ad- 
ditional powers of^ In T/g^ ~ In /3 which might 
also build up to a power of /3. Such an ambigu- 
ity is of the same order as the deviation of the 
1//3 extrapolation from the leading perturbative 
result in Fig. 3. We checked this extrapolating 
with l/l3^-', < e < 0.5.4 



■^Such logarithms are due to the IR divergencies specific to 
d=2 + l and appear, e.g., in the calculation of the Debye 
mass in this theory [10]. 

^This results in smaller extrapolated ao values. 
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Figure 3. Comparison of the lattice results from 
the linear 1//3 extrapolation to /3 = cxd at fixed Lt 
and from perturbation theory for ao in (5), (6). 

In conclusion, while our work is at a prelimi- 
nary stage both analytically and numerically, we 
can already see that perturbation theory does in- 
deed describe the properties of the interface very 
accurately. 
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